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TOPOLOGY-PRESERVING DIFFUSION OF DIVERGENCE-FREE 
VECTOR FIELDS AND MAGNETIC RELAXATION 
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$_, ■ Abstract 

<f | The usual heat equation is not suitable to preserve the topology of divergence-free 

vector fields, because it destroys their integral line structure. On the contrary, in the fluid 
mechanics literature, on can find examples of topology-preserving diffusion equations for 
divergence-free vector fields. They are very degenerate since they admit all stationary 
solutions to the Euler equations of incompressible fluids as equilibrium points. For them, 

<]^ we provide a suitable concept of " dissipative solutions" , which shares common features 

with both P.-L. Lions' dissipative solutions to the Euler equations and the concept of 
"curves of maximal slopes" , a la De Giorgi, recently used to study the scalar heat equation 
in very general metric spaces. We show that the initial value problem admits such global 
solutions, at least in the two space variable case, and they are unique whenever they are 
smooth. 

> 

(N ■ 1. Introduction 

Related to the numerous literature devoted to "topological fluid mechanics" (see [U 
[TUl [T7J [Hi [TH [20] and many others), there are some highly non-linear (and degenerate) 
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diffusion equations for divergence-free vector fields. A typical example is 

(1.1) d t B + V • (B <g> v - v <g> B) = 0, v = WV-(B® B), V • B = 0, 

where P denotes the L 2 projection onto divergence- free vector fields. Following [T7], we call 
them "magnetic relaxation equations" (MRE). They have a (somewhat artificial) physical 
interpretation as they (tentatively) describe a friction dominated model of incompress- 
ible magnetohydrodynamics (MHD) (in some lose sense, "MHD in porous media"). But, 
they have a specific mathematical interest because of their link with the Euler equations 
of incompressible fluids and the theory of "topological hydrodynamics", as discussed now. 

Let us consider the MRE f 1 1.1 1) on the flat torus D = M. d /Z* d , just for simplicity, and 
sketch their three main properties. First, these equations admits an interesting dissipa- 
tion property for the "magnetic energy" , namely 

j /■ | d|2 r 

-dx+ / |PV- {B® B)\ 2 dx = 



dt J 2 

(this is, of course, formally obtained by elementary calculations, assuming B to be 

smooth). Next, the "equilibrium states", for which the energy no longer dissipates, are 

precisely all stationary solutions to the Euler equations of homogeneous incompressible 

fluids, namely 

PV • (B ® B) = 0, V-B = 0. 
l 



Last, equations (jl.ip are "topology-preserving" in the sense that there is a velocity field 
v = v(t,x) £ M. d that "transports" B, i.e. 

(1.2) d t B + V -(B®v-v®B) = 0. 

What we mean by "topology-preserving" is that, for each fixed time t, the "magnetic 
lines" of B (i.e. the integral lines s — > £(s) satisfying £'(s) = i?(£, £(s)), where t is 
"frozen") are "transported" by the flow associated to v , as time evolves. Thus, these lines 
keep their topology unchanged during the evolution, in particular their knot structure. 
[More precisely, let us use "material coordinates" (t,a) — > X(t,a), so that 

d t X(t, a) = v(t, X(t, a)), X(0, a) = a. 

Then, under suitable smoothness assumptions on v and B, the transport equation (11.2J) 
exactly means Bi(t, X(t, a)) = ^ ■ djXi(t, a)Bj(0, a), i — 1, • • -, d (to check the formula, 
just differentiate it in t and use the chain rule.) This implies that the magnetic lines of 
B at time t are the images by X(t, ■) of those of B at times 0. (To check this statement, 
differentiate X(t,£(s)) in s for each integral line £ of B at time 0.)] 
To summarize these three properties, we can say that the "magnetic relaxation equations" 
(II. ip are a good candidate to drive, as time goes to infinity, each given initial magnetic 
field of prescribed topology to a stationary solution of the Euler equations with the same 
topology. This is clearly part of the more ambitious program of "topological hydrody- 
namics", as developed in the papers of Moffatt [17] and the book of Arnold and Khesin jl]. 

Let us emphasize that the standard linear diffusion equation (on the flat-torus) for 
divergence-free vector fields reads d t B = AB and is certainly not "topology-preserving" 
since it cannot be written as a transport equation (11.21) for any vector field v = v(t,x). 
This is in sharp contrast with the standard linear heat equation for positive density fields 
d t p = Ap, which can be easily put in "transport" form 

(1.3) dtp + V • (pv) = 0, v = -V(\ogp). 

Recently, the heat equation has been studied by Gigli, Gigli-Kuwada-Ohta, Ambrosio- 
Gigli-Savare [HI [TJJ [2], in a very general class of metric spaces. Their method is based on 
the following very simple and remarkable idea (that combines the concept of "curves of 
maximal slopes" introduced by the De Giorgi school for "gradient flows" and the interpre- 
tation by Kinderlehrer, Jordan and Otto [T3], [1] of the heat equation as the gradient flow 
of Boltzmann's entropy for a suitable Monge-Kantorovich metric on the set of probability 
measures). We first say that a pair of measures (p(t,x) > 0,q(t,x) £ M d ) is admissible if 
it solves the "continuity equation" 

(1.4) d t p + V ■ q = 0. 
Next, we formally get, for each admissible pair 
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This suggests to characterize the solutions of the heat equation precisely as those admis- 
sible pairs (p, q) that achieve inequality 

(i.5) 2 |/ plogp+ /M! + |^!< . 

This very simple formulation is quite powerful. First, the set of solutions (p, q) for a 
given initial condition is convex and weakly compact. (Notice that the three function- 
al involved in ( 11. 5p are all convex in (p,q), with possible value +00, the first one being 
strictly convex.) Next, uniqueness of solutions directly follows from the strict convexity of 
p — > p log p. (Indeed, the average of two distinct solutions would lead to a strict inequality 
in (II. 5p . which turns out to be not possible). Finally, formulation (ll.5p makes sense in a 
very large class of metric spaces [HI [121 I2J- Notice that this strategy can also be carried 
out for a rather large class of non- linear diffusion equations, as explained in Appendix 1. 

We have spent several lines in explaining this recent approach to the scalar heat equa- 
tion precisely because we are going to follow a similar (but much less successful) way to 
address the (more) challenging magnetic relaxation equations (11.11) . First, we call admis- 
sible solution any pair of time-dependent divergence-free vector fields (B,v) solving the 
transport equation (ll.2p . namely 

(1.6) d t B + V-{B®v-v®B) = 

(which, unfortunately, is non-linear but, at least, has a nice "div-curl" structure in the 
two-dimensional case d = 2). Next, we observe that, for any (smooth) admissible pair 
(B,v): 

d 

di 
(just because of (II. 2p ) 



\B\ 2 dx = 2 I v V- (B®B)dx. 
= 2 / v ■ PV • (B ® B)dx 



(because v is divergence-free) 

= j \v - PV • (B <g) B)\ 2 dx - I \v\ 2 dx - I |PV • (B <g> B)\ 2 dx. 

Therefore, we can characterize the solutions (B,v) of (11. ip . just by asking them to be 
admissible and satisfy the following inequality 

— / \B\ 2 dx + / \v\ 2 dx + / |PV • (B <g> B)\ 2 dx < 0. 

Unfortunately, with respect to the simpler case of the scalar heat equation, we loose 
on two sides. First, the transport "constraint" (II. 6p is not linear (in contrast with the 
continuity equation (11.41) ). Second, the energy inequality involves a non-convex functional 
of B, namely J |PV ■ (B ® B)\ 2 dx. However, in the present paper, we will be able to 
overcome some of these difficulties. The output is a "reasonable" concept of "dissipative 
solutions" , sharing the same strength and weakness as Lions' dissipative solutions to the 
Euler equations [H]: weak compactness (at least in the two-dimensional case d = 2, 
in our case) and uniqueness whenever they are smooth. To finish this introduction, let 
us mention the analysis by Nishiyama of related magnetic relaxation equations [HI [19], 



based on the concept of measure- valued solutions, as well as two recent papers [8|, [9] on the 
concept of dissipative measure- valued solutions for elastodynamics and fluid mechanics. 



2. Dissipative solutions to the magnetic relaxation equations 

Preliminaries. Let D be the d- dimensional flat torus (IR/Z) d , and 1 1 • 1 1 be the L 2 -norm on 
D. Let us denote L^ iv (D; M d ) the space of all L 2 vector fields on D that are divergence- 
free 
Next we define 

(2.7) L{B) = \\FV ■ (B ® B)\\ 2 e [0,+oo], 

for all B G L\ vQ (D\ M. d ), where P denotes the (Helmholtz-Leray) orthogonal projector 
L 2 — >L 2 diY(i (D\ M. d ). A more precise definition is obtained by duality 

L(B) = sup (B®B): (Vz + Vz T )dx - \\z\\ 2 , 

(2.8) Jd 

zeC l (D- R d ), V -z = 0, 
where (B <g) B) : (Vz + V ' z T ) should be understood as 



d 

J2 BiBj&Zj + a 



■jZi). 



L is not a convex function of B and we have to find a substitute for L. This is why, for 
each real nonnegative number r e M + , we define, for B e L^ iv (D; R d ) 



K r {B) = sup{ / (B®B): [Vz + Vz T + rl)dx - \\z\ l 
Jd 



l\ ,■[ L) I ■ hllDi 7 I U ' ' L) I . \ \ „ -1- \ .:. I I III .1 - I I ... I 

(2.9) 

2 G C'(D ; M" 1 ), V-^ = 0, V^ + Vz T + rJ>0} 



(in the sense of symmetric matrices, where / denotes the identity matrix). Notice that 
K r (B) is always bounded from below by r||5|| 2 (take z = in its definition) and is a 
convex function of B (as a supremum of positively curved quadratic functions of B, thanks 
to the constraint Vz + Vz T + rl > 0). In addition, we can recover L out of the K r since 



sup K r (B) - r\\B\\ 2 = sup{ / (B ® B) : (Vz + Vz T + rl)dx - \\z\\ 2 - r\\B\\ 2 , 

r>0 Jd 

r > 0, z <E C^D ; R d ), V • z = 0, V^ + Vz T + rl > 0} 

(2,10) = sup{ I (B®B): (Vz + Vz T )dx - \\z\\ 2 

Jd 

r>0, zeC l (D;R d ), V ■ z = 0, Vz + Vz T + rl > 0} 

= L(B). 



5 
Definition of dissipative solutions to the magnetic relaxation equations. 
Definition 2.1. Given a final time T > and B° e L^ iv Q (D; M. d ), we say that a pair 

(t,x) E [0,T] x D->(B,v){t,x) ER d xR d 

is a dissipative solution of the MRE U.l\) with initial condition B° , if 
i) B is weakly continuous from [0, T] to L"^ iv0 (D; IR d ) with 5(0) = B° ; 
ii) v is square-integrable from [0, T] to Lj iv0 (D;R d ); 
Hi) B, v solves the transport equation Al.fy) , namely 

d t B + V • (B <g> v - v <g) 5) = 0, 

in i/ie sense o/ distributions; 
iv) 

\\B(t,.)\\ 2 + Jj\\v(s,-)\\ 2 + K r{s) (B(s,-))}exp(R(t) - R(s))ds 

<||B(0,0|| 2 exp(i2(t)) VtG[0,T], 
/or a// nonnegative function t — >■ r(t) > 0, wii/i -R(t) = /„ r(s)ds, where K r is defined by 

3. Stability of smooth solutions among dissipative solutions 

We ignore whether or not the MRE (11.11) are locally well-posed in any space of smooth 
functions. (As a matter of fact, this is a very interesting issue, since these equations can 
be considered as parabolic only in a weak sense.) However, we can prove: 

Theorem 3.1. Assume D = (IR/Z) d . Let (B,v) and (J3,w) be respectively a dissipative 
and a smooth solution to the MRE U.l\) up to time T. Then, there is a constant C 
depending only on the spatial Lipschitz constant of (/3, uj), so that, for all t G [0, T], 

(3.12) 



\\(B-P)(t,-W + J exp(C(t-s))-\\v( S ,-)-co(s,-)\\'ds <\\(B-(3)(0,-)\\'exp(Ct) 

This implies the uniqueness of smooth solutions among all dissipative solutions, for any 
given prescribed smooth initial condition. 

Proof. For simplicity, we use notations B t , v t for B(t,-), v(t,-), etc... Let us introduce 
for each t E [0, T] 

iV t = ||Bo|| 2 exp(rt)- f (\\v s \\ 2 + K r (B 8 ))exp(r(t - s))ds 

Jo 

where r is a nonnegative constant to be chosen later. By definition, we have 

(j t -r)N t + \\v t \\ 2 + K r (B t ) = 

(in the distributional sense and also for a.e. t e [0,T]). Since (B,v) is a dissipative 
solution we get from ( 11.21) 

N t >\\B t \\ 2 ,Vte[0,T\. 
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We now set 

e(t) = \\Bt- A|| 2 + (N t - \\B t \\ 2 ) = N t - 2{{B t ,0t)) + ||A|| 2 , Vt G [0,T]. 

where ((•, •)) denotes the L 2 inner product, and compute the time derivative e'(t) of e(t). 
We already know that 



d 



N t = rN t -\\v t \\ 2 -K r (B 



dt 
Next, since (B,v) is a dissipative solution we get from ( 11.2ft 

-((s t , A)) = / s«C8«, t + %(A« - /^)) 

where we use notations f3 U j = dj((3 t )i, etc. ..and skip summations on repeated indices i, 
j... Thus, we find 

e'(t) = rN t - \\v t \\ 2 - K r {B t ) + J 2{j3 - B) u f3 u ,t - 2B t{ v tj {j3 titj - f3 tj>i ). 

By definition ( 12.9ft . we can find a constant r, depending only on the spatial Lipschitz 
constant of to, large enough so that (by setting z = —u in (12. 9ft ) 



-K r (B t ) < / B ti B tj (ujtij + uj t j,i)dx -r\\B-, 



, ,2 , II 1 12 



Thus 

e'(£) < r(iV t - ||B t || 2 ) - H^H 2 + m|| 2 + J t 
where 

Jt = / B ti B tj (uJtij + Wtj^) + 2(/3 — B) ti ft ti)t — 2B ti v t j{/3u,j — Ptj,i)- 

We may write 

where J t , J 4 L1 , J 4 L2 , J t c are respectively quadratic, linear, linear, and constant with respect 
to B — f3 and v — u, with coefficient depending only on uj, j3: 

J? = f(B - f3) ti (B - 0) tj (u ti ,j - wtj-i) - 2(B - /3) ti (v - u) tj {P*j ~ tj ,i) 



J t L1 ----- I 2[B - f3) ti [(3 t] {ujtij + u tj ,i) ~ Pu,t - u)tj(pHj ~ Ptj,i)] 
J t L2 = [2(v-u) tj /3 u /3 tj , 



Let us reorganize these four terms. By integration by part of its second term, we see that 
J t c = 0, using that (3 and uj are divergence-free. Next, since B t — f3 t is divergence-free, we 
have 

J t L1 = J 2(B- P)ti[-Pu,t - WtjPtij + PtjWtij] 
and we may reorganize 

J t L2 = 2((v t - u t , ut)) + 2(v- w)tj\ptiptj,i ~ u t j] 



Thus 

e'it) < r(N t - \\B t \\ 2 ) - \\v t - io t \\ 2 + J? + J t L 
where 

J 4 L = -2((B t - fr, dtPt + (u t ■ V)A - (A • V)w*)) + 2((v t - u t , V(A ® A) - wt)). 
Clearly 

for some constant C depending only on the spatial Lipschitz constant of ((3,u). Finally, 
we have obtained 

e\t) + ±\\vt - cu t \\ 2 < r(N t - \\B t \\ 2 ) + C'\\B t - f3 t \\ 2 + J t L , 

and, therefore, by definition of e(t) 

e(t) = \\B t -/3 t \\ 2 + (N t -\\B t \\ 2 ) } 
we get 

e'(t) + -\\v t -oo t \\ 2 <Ce(t) + J t L , 

where C is another constant depending only on the spatial Lipschitz constants of (/3, u) 
(through R). By integration we deduce 



e(t) + 



f e ^ c -\\v s - co s \\ 2 ds < e(0)e tc + f e^ c J^ds. 
Jo 2 J 



Finally, let us remember that e(t) > \\B t — /3 t || 2 with equality for t = (since N t > ||-Bt|| 2 
with equality at t — 0). Thus, we have shown 

Lemma 3.2. Assume D = (IR/Z) 01 . Let us fix T > 0. Let (B,v) be a dissipative solution 
of the MRE U.l\) up to time T , and ((3, uj) be any pair of smooth functions ((3, to) chosen in 
-^div o{D] B^)- Then there is a constant C depending only on the spatial Lipschitz constant 
of (J3,u), up to time T, so that, for all t e [0, T], 

(3-13) 

\\B t -/3 t \\ 2 + f t e (t - s)C l\\v s -u s \\ 2 ds<\\B -/3 \\ 2 e tC + /*V^J>, 
Jo * Jo 

J 4 L = -2((B t - (3 t , dtfit + (u t ■ V)A - (A • V)w t )) + 2((v t - u t , V(f3 t ® ft) - u t )). 
We now see that J 4 L exactly vanishes as {(3,u>) is a smooth solution to the MRE fll.ip . 
dtP + (u-V)P-(P-V)u = 0, w = PVC9®/3), V-£ = V-w = 0, 
and the proof of Theorem 13.11 immediately follows. 



4. Existence of dissipative solutions in two space dimensions 

Theorem 4.1. Assume d = 2 and D = (R/Z) 2 . Let T > and fix an initial condition 
B° e L diy (-D; M d ). T/ien i/iere is at least one dissipative solution (B, v) of the MRE U.l\) 
up to time T. This solution can be obtained as the limit in C°([0,T]; L div0 (D; M. d ) w ) x 
L 2 ([0,T]; Lj iv0 (D; M. d )) w , as parameters (e,/i, v) go to zero, of the unique solution of the 
MHD system (with friction and viscosity) 

e(d t v + V-v®v) + v + Vp- uAv = V • B®B, 
(4.14) 
v ; d t B + V- {B®v-v® B)-uAB = 0, V • v = V ■ B = 0, 

with smooth initial conditions chosen so that B(0, ■), y/e v(0, •) approach respectively B° 
and in L 2 . 

Proof. Since d = 2, the MHD system has global smooth solutions for any smooth intial 
condition [21] . From these MHD equations, respectively multiplied by v and B and 
integrated over D, we get two straightforward estimates 



(415) i iw 



/ -ff + / M 2 + ^ / I Vv \ 2 = ~ f(B®B): Vv 

+ u [ \VB\ 2 = f{B®B) :Vv 



dt J 2 
We first add up these estimates in order to get the total energy balance 

(4.16) 4- I |i?|2 + eH2 + / \v\ 2 + H [ I Vv| 2 + v [ \VB\ 2 = 0. 



dt J 2 

This implies that B and v are respectively uniformly bounded (with respect to (e, /i, v)) 
in L°°([0,T\,L 2 (D)) and L 2 ([0,T], L 2 (D)). Using the second equation of gUD, we also 
see that B is uniformly bounded in C 1/2 ([0,T], (Lip(£>; JR d )') (where Lip(D;R d )' denotes 
the dual of the space of vector- valued Lipschitz functions). Indeed: 

(4.17) 

| I dt' f[(B(t', x) <g) v{t', x)) : (Vz T - Vz){x) - uVB(1/, x) : V ' z{x)}dx\ 2 
< \t-s\ Lip(z) 2 [4 sup f \B{t'\x)\ 2 dx f \v{t\x)\ 2 dxdt' + v 2 f \VB{t\x)\ 2 dxdt'}. 

t"6[0,T]7 J J 

We deduce that, as e,/j,,u go to zero, B and v are compact respectively in the spaces 
C°([0,T},L 2 divfi (D;R d ) w ) and L 2 ([0,T], L 2 divfi (D;R d )) w , where subscript w refers to the 
weak topology of L 2 . 
Next, we use the weak formulation (in x) of the first equation of (I4.14I) . namely: 

e(d t v + V-v®v)+v + Vp- fiAv = V ■ B®B, 

to get, for any fixed smooth test function z(t, •) valued in L div (D; M. d ) 



Adding up the energy balance (14.161) . namely: 



d fAB\ 2 + e\v\\ 



we obtain 



dt.l 2 

d [,\B\ 2 + e\v\ 2 



+ / \v\ 2 + /2 / \Vv\ 2 + u / |VB| 2 = 0, 



(4.18) 



,, , . ,, + eviZi) - I (e(viZ ijt + ViVjZ i:j ) + /jz id v id ) 

/i |2 112 r I I 12 

(H- + S^i - H.) + y (i^J- + A*l V^| 2 + H V5| 2 ) = 0. 

Let us now introduce any nonnegative function t — >■ r(t) and i?(t) = J r(s)ds, and assume 
that (,2jj +%-,i + r(t)<5jj) is a positive matrix. After multiplication by 2, we may rearrange 
( 14~T8|) as 

(j t ~ r ) (\ B \ 2 + e \ v + z \ 2 ) + / (M 2 + ^ 5 j(^j + Hi + r6 ij) - 

(4.19) +r / e|f + z| 2 - — / e\z\ 2 -2 (e(viZu + ViVjZij) + +/j,z i> jV i>j ) 



\z\ 2 ) 



\v + z\ 2 + 2fi\Vv\ 2 + 2z/|V5| 2 ) = 0. 



Thus, 



(4.20) (j t - r) f(\B\ 2 + e\v + z\ 2 ) + j {\v\ 2 + BiBfaj + z hi + r5 tj ) - \z\ 2 ) < V (t) 



where rj(t) depends on the fixed test function z and goes to zero in I/ 1 ([0, T]) with (e, /i, is), 
since v is uniformly bounded in L 2 . Next, we integrate in time this differential inequality, 
and, then, we let (e, /i, v) go to zero, assuming that the initial condition B(0,-) and 
\/ey(0, ■) converge in L 2 respectively to the given initial condition B° and to 0. After 
these operations, we obtain for any accumulation point of the (B,v), still denoted by 
(B,v), 

\\B(t,-)\\ 2 + [ \\v(s r )\\ 2 ex P (R(t)-R(s))ds 
Jo 



(4.21) 



/ / dx{B i Bj(2zi i j + rSij) — \z\ 2 )(s,x) exp(i?(t) — R(s))ds 
Jo Jd 

<\\B(0,.)\\ 2 exp(R(t)) Vte[0,T], 

using the positivity of (zij + Zj ti + rS^). Next, taking the supremum with respect to z, for 
fixed r, leads to the dissipation inequality ( 12. lip involved in the dissipative formulation of 
the magnetic relaxation equations. However, we are still left with the problem of passing 
to the limit in the transport equation ( 11.21) . We do not see any clue for that, except in the 
bidimensional case d = 2, where we have a nice "div-curl" structure. Indeed, as d = 2, at 
least locally, we can write 

B = {d 2 A,-d l A), 

for some scalar potential A{t,X\,Xo) G M. Then, the transport equation (11.21) can be 
integrated out as 

d t A + V ■ (Av) = 0, 
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and we can pass to the limit, since |VA| = \B\ and v are well controled in L 2 (using 
estimate (I4.17P to handle the time dependence). This concludes the proof. 

5. Appendix 

5.1. A general framework for dissipative equations. A rather general framework 
that one can encounter in several situations of Mechanics and Physics is as follows. (We 
do not claim any novelty in it, see, for instance, [E] for somewhat related issues.) Working 
on the d- dimensional flat torus D = M. d /Z d , for simplicity, we call admissible a pair (B, E) 
made of two time-dependent differential forms, one of degree k, say B, that we assume to 
be closed dB = 0, and one of degree k — 1, say E, linked by 

(5.22) d t B + dE = 0. 

(A simple example being k = d, as in the scalar heat equation where B = p, E = q.) Next, 
we are given a scalar function L(E, B), called "Lagrangian" , that we suppose convex in E 
(with possible value +oo). Then, we define its Legendre-Fenchel transform with respect 
to E, called "Hamiltonian' 

(5.23) H(D,B) = supE-D-L(E,B), 

E 

and introduce the "defect" function 

def (B, E, D) = L{E, B) + H(D, B)-E-D>0 

(5.24) with equality if and only if E = dxH^D, B). 

if and only if D = d 1 L(E,B). 

Now, we are given a convex function 9 and compute for any (smooth) admissible pair 
(E,B) 

— 1 6(B)dx = - ! 6\B) -dE = - I E-D 

where D = 5(6'(B)) (here 5 = (— l) fe_1 * _1 d* is the Hodge co-differential). Thus 

— J 6{B)dx = /(def (5, E, D) - L{E, B) - H(D, B)) > - f{L{E, B) + H{D, B)) 

with equality if and only if E — diH(D,B) (pointwise). This suggests a "dissipative 
formulation" for the (highly) non-linear equation 

d t B + dE = 0, 

(5 ' 25) E = d 1 H(5(6'(B)),B) 

We call dissipative solutions of this equation any admissible pair (E, B) such that 

d 



(r) ., 0) ijf , 0(B)dx+ l(L(E,B) + H(D,B))<0, 

where D = 6(0'(B)). 
Of course, if 

r (L(E,B)+H(5(9'(B)),B)) 

turns out to be a convex function of the pair (E, B) the analysis gets rather simple, but 
there is little chance to get interesting examples of this type unless k = d. As a matter 
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of fact, in the case k = d, denoting (E,B) = (q,p), we get the rather general non-linear 
scalar diffusion equation 

d tP = V-d 1 H{V(6 , (p)),p). 

In the special case L(q,p) = pc(-), we get H(v,p) = pc*(v), where c is a convex function 
and c* its Legendre-Fenchel transform. The resulting equation reads 

d t p = V-( P Vc*(V(9\p))). 

(Notice that this equation can be also handled by optimal transport methods [22], using 
"cost function" (x,y) — > c(x — y) and "entropy function" 9.) The further choice 

9(p) = plogp, c(w) = \w\ 2 /2 

leads to the linear heat equation. Another example is the "relativistic heat equation" 
[31 [15] , for which 

9(p) = plogp, c(w) = -a/1 - \w\ 2 . 

5.2. A "topology-preserving" diffusion equation for divergence-free vector fields, 
based on Born-Infeld Electromagnetism. In order to provide a non-scalar applica- 
tion of the general framework, let us consider the special case d = 3 and k = 2. So B 
is a closed 2— form in three space dimensions, which corresponds to a divergence-free 
vector field, while E is a 1— form, i.e. a vector field. We use classical notations x for the 
wedge product as well as Vx for d, which here is the curl operator. So, we call admissible 
solutions any pair of fields (E, B)(t, x) G M 3 satisfying 

(5.27) d t B + V x E = 0. 

Let us introduce the Born-Infeld "Lagrangian" [5j, parameterized by constant A > 0: 

(5.28) L X (E, B) = -a/A 2 + \B\ 2 - \E\ 2 - \- 2 (E ■ Bf. 

Function L\, for each fixed value of B is convex in E (with infinite value as the term 
under the square root gets negative). The Legendre-Fenchel transform with respect to E 
can be easily computed and is given by the " Hamiltonian" 



H X (D,B) = J\ 2 + \B\ 2 + \ 2 \D\ 2 + \Dx B\ 2 
(5 29) Z_ 

= y/{X 2 +\B\ 2 )(l + \D\ 2 )-{D-B) 2 . 

Let us now consider a convex function 9 : R 3 — > R. From the general framework, we get 
a "dissipative formulation" for the (highly) non-linear equation 

(5.30) d t B + VxE = 0, E = d 1 H x (Vx(6'(B)),B) 

by calling dissipative solution any admissible pair (E, B) such that 

j t J 9(B)dx + f(L x (E, B) + H X (D, B)) < 0. 

( 5 - 31 ) where D = V x (9'(B)) 

H X (D,B) = x/(X 2 + \B\ 2 )(1 + \D\ 2 ) - (D ■ B) 2 . 
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The Born-Infeld Lagrangian has two remarkable properties [5J [7]: as A — x oo, we recover 
the classical Maxwell Lagrangian for electromagnetism (as the limit of (A + L\)\)\ as 
A —)• 0, we get 
(5.32) 

H (D, B) = y/\B\\l + |£>| 2 ) - (D ■ B) 2 , L (E,B) = -y/\B\ 2 - \E\ 2 , E ■ B = 0, 

which, interestingly enough, includes the pointwise constraint E ■ B = 0. This exactly 
means there is a vector v G M 3 such that E = B x v and v can be defined, for instance, 
by setting 

Ex B 

v = 1bT- 

So, the constraint dtB + V x E = becomes 

at5 + Vx(Bxti) = 0, 

or, equivalently, 

d t B + V • {B ® v - v ® B) = 0, 

which is exactly the "topology-preserving" equation ( II .2p (with, here, a vector field v that 
is a priori not divergence- free). Let us now express v in terms of B. We have 

x D\B\ 2 -(D-B)B 
E = dl H (D,B) = L_L__ 

where D = V x (^'(5)). Thus 



|S| 2 ffo(AB) y/\B\* + \DxB\ 2 ' 
So, in the case A = 0, equation (I5.3QJ) can be rephrased as 

Sifi + V • (B ® v - u ® B) = 0, 

DxB 

(M3) " = ^|B|» + |JxB|» 

D = Vx (0'(B)). 

The resulting system looks very much like the magnetic relaxation equation (11.11) (without 
divergence-free constraint for v). For instance, in the case 9(B) = \B\, we get 

B®B 



V- 



Ho 



where H = ^\B\ 2 + \D x B\ 2 and D = V x (-B/|S|). 
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